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Abstract We revisit the vertical structure of neutrino- 
dominated accretion flows (NDAFs) in spherical coordi- 
nates with a new boundary condition based on the me- 
chanical equilibrium. The solutions show that NDAF is 
significantly thick. The Bernoulli parameter and neu- 
trino trapping are determined by the mass accretion 
rate and the viscosity parameter. According to the dis- 
tribution of the Bernoulli parameter, the possible out- 
flow may appear in the outer region of the disk. The 
neutrino trapping can essentially affect the neutrino ra- 
diation luminosity. The vertical structure of NDAF is 
like a "sandwich" , and the multilayer accretion may ac- 
count for the flares in gamma-ray bursts. 
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1 Introduction 



black hole 



Neutrino-dominated accretion flows (NDAFs) involve a 
stellar black hole of 2 ~ 10 Mq, accreting wih a hy- 
percritical rate, in the range of 0.01 ^ 10 Mq s~^. In 
the past decade, a number of st u dies h a ve investigated 
this m o del (e.g. Popham et al. J 19991): Naravan et al 



2001) ; |Kohri fc Mmeshige v^^^m, i-i^i iv^attcu cu solution was establislied. llius. \Lm et al.l puiUal l pre- 

2002) ;|Rosswog et al.| (|2003|);|Kohri et al.| (|2005|):|Lee et al.|[;nted the vertical struct ure of NDAF calc ulating with 



(I2OO2I): bi Matteo et al 



In cylindrical coordinates. iGu &: Lul (|2007[ ) discussed 
the underlying importance of taking the explicit form 
of the gravitational potential for calculating slim disk 
solutions, and pointed out that the Hoshi form of the 
potential ( Hoshil Il977l ) is vali d only for geometrically 
thin disks with H/R < 0.2. IlIu et al.l pOOSi ) found 
that NDAFs have both a maximal and a minimal pos- 
sible mass accretion rate at their each radius, and pre- 
sented a unified description of all the three known 
classes of optically thick accretion disks around black 
holes, namely Shakura-Sunyaev disks, slim disks, and 
NDAFs. These two works are, however, based on the 
simple on e-zone ver t ical h ydrostatic equilibrium. Fur- 
thermore, IGu et al.l (|2009l ) revisited the vertical struc- 
ture of advection-dominated accretion disks in spheri- 
cal coordinates and showed that those disks should be 
geometric ally thick r ather than being slim. As a con- 
tinuation, ICai et all ( 2010t ) found an analytic relation 
Cbo/vk& = [(7 — 1)727]^^^ for the vertical mechanical 
equilibrium, where Cso is the sound speed on the equa- 
torial plane, vk is the Keplerian velocity, O is the half- 
opening angle of the fiow, and 7 is the adiabatic index. 
However, the detailed radiation cooling was not consid- 
ered in that work, and therefore no therma l equilib rium 



olution was established. Thus, iLiu et al.l (|2010al ) pre 



( 20051): IGu et a l.' (200 ^: IChen fc Beloborodovl (IzOOTI) : 
Janiuk et al.rTsOO?, '): .Liu et al.l (|2007i boOSl boiOallbh V 
The model can provide a good understanding both 
the energetics of gamma-ray bursts (GRBs) and the 
processes of making the relativistic and baryon-poor 
fireballs by neutrino annihi lation or magn e tohyd rody- 
namic processes (see e.g. IPopham et al. (1999) and 
Di Matteo et al.l (|2002l ) for references). 
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the numerical method of Gu et al. ( 20091 ). We found 
the luminosity of neutrino annihilation is enhanced by 
one or two orders of magnitude. The empty funnel 
along the rotation axis can naturally explain the neu- 
trino annihilable ejection. 

The Bernoulli parameter, defined as the summa- 
tion of kinetic energy, potential energy and enthalpy 
of the accreting gas, is an essential quantity in accre- 
tion flows since it can be used to measure the possibil- 
ity o f arising outflows or winds ( Naravan fc Yil 1994 . 
I995I) . It is well-known that the positive Bernoulli 



parameter is necessary for the occurrence of outflow. 
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Naravan fc Yil (119951) considered rotating spherical ac- 
cretion flows ranging from the equatorial plane to the 
rotation axis. They found that the Bernoulli param- 
eter is positive in advection-dominated conditions, es- 
pecially close to the rotating axis. Here we also can 
evaluate the vertical distributions of Bernoulli param- 
eter basing on our solutions, and give a judgement for 
whether possibly arising outflows. 

Moreover, the neutrino trapping is a nother no- 
table point in the studies of NDAFs. iBegelman 



2 Equations and Boundary Conditions 

2.1 Equations 

We consider a steady state axisymmetric accretion flow 
in spherical coordinates (r, 9, 0), i.e., d/dt = d/d4> = 0. 
We adopt the Newtonian potential -0 = ~GM/r, since 
it is convenient for self-similar assumptions, where M 
is the mass of the central black hole. The basic equa- 
tions are composed of c ontinuity, monientum and en 



( 19781) m-esen t ed photon trapTj ing i n Bondi accretion. ergy e q uations, {see e.g. Xue fc Wanj (|2005|) ; |Gu et al 



Ohsuea et all (|20Q2l , l20Q3il20CI5[ ) and lOhsuga fc Mineshigc| (|2009l) : lU u et all (l2010al)) 



( 201lh discussed the photon tr apping in black h ole ac- 
creting systems. For NDAFs, lOi Matteo et al.l (j2002) 
argued that the neutrino trapping occurs when the neu- 
trino optical depth is larger than 2/3, and the inner 
regions of their solutions are all opaque for neutrinos 
even for the moderate accretion rate due to the using 
of Newtonian potential. Thus, in their solutions, neu- 
trinos are significantly trapped by the fiow, and the 
neutrino ann ihilation lumino sity becomes too low to 
boost GRBs. IGu et al.l (|2006l ) found that the neutrino 



optical depth would significantly decrease and there 
would be no any neutrino trapping if the general rela- 
tivistic effect is considered. Accordingly, they argued 
that NDAFs still have eno ugh luminosity to work as the 
central engine of GRBs. IChen fc Beloborodovl (|2007l) 
discussed the neutrino trapping radius in NDAFs. They 
found that this radius depends on the accretion rate, 
viscosity parameter and spin of black hole. These works 
all concentrated on discussing neutrino trapping in the 
equatorial plane of disks, so it is necessary to study the 
neutrino trapping in the vertical structure of NDAFs. 

Our purposes of this paper are to investigate the 
vertical structure, distributions of Bernoulli parameter 
and neutrino trapping in NDAFs with detailed neutrino 
radiation and new boundary condition. In section 2, 
we show the numerical results obtained by solving the 
vertical differential equations with the self-similar as- 
sumptions in the radial direction and a boundary con- 
dition on disk surface, then we give the definition of 
Bernoulli parameter and the criterion of neutrino trap- 
ping. In section 3, we show the vertical distributions of 
the Bernoulli parameter, neutrino trapping and other 
physical quantities, such as radial velocity and electron 
fraction of different cases. The possible applications in 
GRBs of these results is also mentioned. Conclusions 
and discussions are made in section 4. 



Similar to iNaravan fc Yil (|l995l ) , we adopt the self- 
similar assumptions in the radial direction to simplify 
the continuity and momentum equations, then reform 
and obtain the following vertical equations: 
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(1) 



(2) 



(3) 



where Vr and are the radial and azimuthal compo- 
nents of the velocity {ve = 0), the sound speed Cg is 
defined as Cg ^ p/ the Keplerian angular velocity is 
Ok = [GM jr'Y^'^ 1 ^^'^ ^ is the constant viscosity pa- 
rameter. 

According the continuity equation, the mass accre- 
tion rate can be written as 



■it/2 



pv, 



■ sin 9d9, 



(4) 



where 9o is the polar angle of the surface and p is the 
density. 

The total pressure of gas is expressed as the summa- 
tion of four pressures: 



P = Pgas + Prad +Pc+P,y, 



(5) 



where Pgas, Prad, Pc, and p^, are the gas pressure from 
nucleons, the radiation pressure of photons, the degen- 
eracy pressure of electrons, and the radiation pressure 
of neutrinos, respectively. Detailed expression s of th e 



pressure components were given in ILiu et al.l (|2007|). 
It should be pointed out that NDAFs are extremely 
optical thick for photons, so the radiation pressure of 
photons at anywhere can be written as 



Prad 



(6) 
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Fig. 1 Variations of the half-opening angles of the disk surface {n/2 — 9o) (solid lines), B — (dashed lines), t — 1 (dotted 
line) with radii r/rg for which the given parameters are set as (Case 1) a = 0.1 and M = 1 Mq s~^, comparing (a); (Case 
2) Q = 0.01 and M = 10 Mq s"S comparing (b). 



where a is the radiation constant and T is the temper- 
ature of gas. Additionally, we assume the polytropic 
relation in the vertical direction, p = Kp^/^ , where K 
is a constant and 7 = 4/3 is the adiabatic index. 

Considering the energetic balance, the energy equa- 
tion is written as 



'adv 



(7) 



where Qvisi Qadv: and Qi, are the viscous heating, ad- 
vection and neutrino cooling rates per unit area, re- 
spectively. In this paper, we ignore the cooling of pho- 
todisintegration of a-particles and other heavier nu- 
clei. The viscous heating rate per unit volume gvis = 
vpr'^[d[Vff,/r)/drY and the advective cooling rate per 
unit volume gadv — pvr(de/dr — [p / p'^)d p / dr) (e is the 
internal energy per unit volume) are expressed, after 
self-similarly simplification, as 



9adv 



9 apvl 



3 (p - Pe)Vr 



(8) 



(9) 



where the entropy of degenerate particles is neglected. 
Thus the vertical integration of Qvis and Qadv are the 
following: 



Qvis = 2 / q^isT sm 

Bo 



(10) 



The coolin g due to the ne utrino radiation can be 
defined as ( Lee et al.ll2005r ) 



= 2 



e 



'r sm 



(12) 



where Qi, is the summation of the cooling rates per unit 
volume due to the Urea processes, electron-positron 
pair annihilation, nucleon-nucleon bremsstrahlung, and 
Plasmon decay (hereafter we represent then with qi 
(i=l, 2, 3, 4), respectively. Different type neutrinos 



involv e different processes, for details, see iLiu et al 
(120071)); = Ta,v + Ts^v IS the ncutrino optical depth 



including absorption optical depth Ta^u and scattering 
optical depth which can be defined as 



E 



3 



ajUjrdO, 



(13) 



(14) 



where Uj and rij {j = 1, 2, 3, 4) are the cross sec- 
tions and the number density of protons, neutrons, a- 
partic le s, and electrons, respectively {e.g. iKohri et al 
( 20051 ): IChen fc Beloborodovl (|2007[ )). In order to em- 
body the complicated microphysics in NDAFs, we need 
to evaluate the electron fraction = ?ip/(np + 71,1) 
by solving equations (42), (44) and (45) in Liu et al. 
(2007). 



Qadv = 2 / gadYT-smf 

'Bo 



(11) 



2.2 Boundary condition 

A boundary condition is required for solving the equa- 
tions. The Eddington luminosity, which is defined by 
the equilibrium between the gravity and the radiation 
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pressure, is closely related to the energy-liberation rate 
associated with mass-accretion processes. If the gradi- 
ent of radiation pressure is larger than the gravity, out- 
flows may occur, so we considered that it could present 
a natural bou ndary condition for the accretion disk. 
Deviated from lLiu et al.l (|2010al ). here boundary condi- 
tion is set to be the mechanical equilibrium, like solving 
principle of Eddington luminosity. Furthermore, we can 
ignore the radial gradient of radial velocity and pressure 
at the surface of the disk if the disk is thick. There are 
three forces to balance, namely gravity, radiation pres- 
sure force and centrifugal force in non-inertial reference 
frame. The mechanical equilibrium can be written as 



Prad \e=Ba ""T 



2GMm,. 



cot(?o, 



(15) 



combined with eq. (6), the surface temperature can be 
defined as 



Tsurf = ^coteio S 



(16) 



where mu is the mean mass of a nucleon, ctt is the 
Thompson scattering cross. 

2.3 Bernoulli parameter and neutrino trapping 

The Bernoulli parameter of th e accr eted matter is ex- 
pressed as {e.g. iNaravan fc Yil \l99^ ) 



3 Numerical Results and Possible Applications 
to GRBs 

Figure 1 shows the variations of the half-opening an- 
gles of the disk surface {tt/2 — Oq) (solid lines), B — 
(dashed lines), i = 1 (dotted line) with radii r/rg 
for which the given parameters are set as (Case 1) 
a = 0.1 and M — 1 Mq s""'^, comparing (a); (Case 
2) a = 0.01 and M 10 Af© s"\ comparing (b). Here 
Tg ~ 2GMj(? is the Schwarzschild radius. In figure 
1(a), the variation of half-opening angle from ~ 0.56 to 
~ 1.25 with r from 3rg to 500rg. The region of B > 
appears at r ~ lO'Tg, and expands continuously with 
the increase of radius. There is no region of neutrino 
trapping. In figure 1(b), the variation of half-opening 
angle from ^ 0.91 to ^ 1.50. The region of _B > 
appears at ^ 9rg, and increases to near the surface of 
the disk at ~ 500rg. The region of neutrino trapping 
occurs from inner area to ~ 46rg. We also notice the 
variation on opening angle of t = 1 is from ^ 0.4 to 
0. We consider that the radius and the open angle of 
5 = and t = 1 is closely detern iined to the accretion 
rate and the viscous parameter (jChen fc Beloborodov 
200i . 

The existence of the B > region implies the possi- 
ble outflow. B < 0, everywh ere is a sufficient conditio n 
for the absence of outflows ( Abramowicz et al. 2000l ). 



B = 



7 2^1. 2^ 2^ GM 
+-{Vr ' 



(17) 



The neutrino trapping can decrease the neutrino radia- 
tion luminosity and annihilation luminosity. In the case 
of Figure 1(b), the neutrino luminosity is obtained as 



which includes the kinetic energy, potential energy and 
enthalpy of accreted gas. 

In classic accretion theory, the radiation energy gen- 
erated near the equatorial pla ne diffuses toward the disk 
surfa ce at the speed of ~ c/3t ( Mihalas fc Weibel Mihala^ 
19841) . where r is the total optical depth. Thus, the 
timescale of radiati ve diffusion is <difi — H/{c/3t) 
(jOhsuga et al.l 120021 ) . where H is the half thickness of 
disk. We get w„ for neutrinos replace c for photons. Ob- 
vious ly, Vn is related t o the e n ergy of neutrino ^ 3.7fcT 
(e.g. iDi Matteo et all (|20oi) : iLiu et al.l |2007l) ). Since 
the accretion timescale iacc is given by —r/vr, the con- 
dition in which the neutrino radiation energy in the 
disk is trapped in the flow and falls onto black hole is 
written as tdiff > ^acc- Here we define a parameter as 



J3r„ 



'■dr. 



(19) 



which decreases 18% in comparison with iLiu et al 
(|2010al) . It is stiU about lO'^'^ergs s^^ in Case 2 to 
produce GRBs according to neutrino annihilation. The 
disk is still thick, thus the volume above the disk shrinks 
and the radiated neutrino density increases. The neu- 
trino annihilation luminosity i s also larger than th e 



results of the p re vious works ( Popham et al. ( 19991) : 



t = 



^diff 
tacc 



rVn 



(18) 



only the optical depth of electron neutrino is calculated. 

We numerically solve equations (1-12) with the 
boundary condition (16) for p, Vr, v^, Cs, p, and T 
by given a, M, M, and r(in this paper, we only con- 
centrate on the stellar black hole with M = 3Mq). 



Pi Matteo et al.l (|2002[ ): iGu et al.l (|2006D : iLiu et al 
( 200?!) ). The 3-dimensional numerical simulations 
have shown that the neutrino annihil ation luminosity 
is like l y to be adequa t e for GRBs (|Ruffert fc Janka 
1 1999f ): Rosswog et al. (2003)). Our model is an ac- 
ceptable estimation in comparison with the simulation 
results. 

Figure 2 shows the variations of the electron fraction 
Yc and the radial velocity Vr with the polar angle 9 for 
which the given parameters are Case 1, comparing (a) 
and (b); Case 2, comparing (c) and (d). The profiles 
of Vr are similar to that of t he optically thin adv ection- 
dominated accretion flows ( Naravan fc Yi|[l995 ). Fig- 
ures 2(a) and 2(c) show the variety of electron fraction 
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Fig. 2 Variations of the electron fraction Yc and the radial velocity Vr with the polar angle 9 for which the given parameters 
are Case 1, comparing (a) and (b); Case 2, comparing (c) and (d). 



in the vertical direction. Electron fraction of Case 1 
near the equatorial plane is larger than Case 2, because 
it is determined by accretion rate more than viscosity 
parameter. Yc is approach 0.5 at the surface of the disk 
in the outer region, which means the matter of the disk 
is almost a-particles and electrons in non-degenerate 
phase. GRBs may be the important contributors to 
local, and possibly galactic, nucleosynthesis, especially 
for rare species such as r-process nuclei. A lot of works 
have studied the nucleosynthesis in the cent ral engine of 
GRBs {e.g. ISurman fc McLaughlinl (|2004f l). The dis- 
tribution of electron fraction in our solutions may be an 
implication for the origin of nucleosynthesis in GRBs. 
In figures 2(b) and 2(d), the radial velocity Vr is obvi- 
ously less than it near the equatorial plane. It has not 
greatly changed between Case 1 and Case 2, because it 
is correlated positively with the accretion rate and vis- 
cosity parameter. The vertical structure o f the disk is 
Hke a "sandwich" {e.g. I Zhang et all (|2000h ). The mat- 
ter near the equatorial plane has been accreted first, 
inevitably, the remnants of the matter near the surface, 
which is almost in non-degenerate phase, has been ac- 
creted later, any radial disturbed motion will destroy 
the remnants to fracture in some parts. Two steps of 
accretion processes with different mass accretion rate 
correspond the GRB and its late flares. We find the 



similar model has been proposed by iLee et al. I (I2OO9I) . 
They found that powerful winds are launched from the 
surface of the disk, driven by the recombination of free 
nucleons into a-particles, which may produce the late- 
time X-ray flares. 



4 Conclusions and Discussion 

In this paper we revisit the vertical structure of NDAFs 
in spherical coordinates and calculate the vertical dis- 
tributions of Bernoulli parameter, neutrino trapping, 
electron fraction and radial velocity. The major points 
we wish to stress are as follows: 

1. We show the vertical structure of NDAFs under the 
mechanical equilibrium, and find that the flow will 
be significantly thick. 

2. The Bernoulli parameter and neutrino trapping are 
determined by the mass accretion rate and the vis- 
cosity parameter. 

3. According to the distribution of the Bernoulli pa- 
rameter, we find that the possible outflow may ap- 
pear, particularly in the outer region of the disk. 

4. The neutrino trapping can essentially affect the neu- 
trino radiation luminosity and the neutrino annihi- 
lation luminosity in turn. 
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5. The distributions of radial velocity and electron frac- 
tion indicate that the vertical structure of NDAF is 
like a "sandwich" , and this multilayer accretion may 
account for the progenitor of the late-time flares in 
GRBs. 

It is generally a two-dimensional problem to study a 
steady and axisymmetric accretion flow, and it is quite 
difficult for directly solving the two-dimensional system 
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(e.a. Pooham et al. (119991): Naravan et al. 


2001 


) 
) 
) 




Kohri & Mineshiee ( 2002 ) : 


Di Matteo et al. 


2002 


Kohri et al. (2005): Lee et al. 


(2005): Gu et al. 


2006 


Chen & Beloborodov (2007): Janiuk et al. (20071 


: Liu et al.l 


(20071)) focused on the radial structure and therefore 



simplified the vertical counter-part of disk by using a 
one-zone approximation. The purpose of the present 
work is, however, to investigate the vertical structure 
in more details. Consequently, it requires us to directly 
solve the differential equations in the vertical direction, 
thus some simplification has to be taken in the radial 
direction. The self-similar assumption may be an ex- 
cellent method to achieve the above purpose. More- 
over, such an assumption has been widely adopted and 
successfully applied to neutrino-radiative disks (e.g. 



Naravan et al. 



(I2OOII)) as well as to photon-radiative 



disks {e.g. iNaravan fc Yil (|1994 )). 

Based on the radial self-similar assumption, we can 
obtain some new res ults compared with the pre vious 
works. For example, IChen fc Beloborodov! (|2007l ) dis- 
cussed the radial distribution of Bernoulli parameter, 
neutrino trapping, and electron fraction, but the verti- 
cal distribution has to be ignored in that work. On the 
contrary, we can present the explicit vertical distribu- 
ti on of physical quanti ties a s shown in the figu res. 



Xue fc Wand (120051) and iJiao fc Wul (120111) studied 



the outflows from accretion disks by assuming V0 ^ Oin 
spherical coordinates. Their results indicated that the 
outflows should be common in various accretion disks 
and might be stronger in slim disks, where both advec- 
tion and radiation pressure are dominant. However, the 
radiation processes are not considered in these works. 
The new vertical structure of NDAFs should be cal- 
culated by reserving vg to discuss the natural outflow 
from the surface of the disk. 

In our next wo r k, we will take neutrino transfer ef- 
fects {e.g. ISawveil (|2003l ): iRossi et al.l |2007l )) into con- 
sideration instead of the polytropic equation of state, 
and calculate the variation of the neutrino luminosity 
and annihilation luminosity with the mass accretion 
rates. The "sandwich" accretion disk model can also 
be studied, and some simulations should be done to 
discuss the possible origin of flares in GRBs. Moreover, 
we will study the complex nucleosynthesis process in 
the vertical direction of NDAFs. 
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